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ABSTRACT 

The MareNostrum Universe is one of the largest cosmological SPH simulation done so far. It consists 
of 1024 3 dark and 1024 3 gas particles in a box of 500 h~ x Mpc on a side. Here we study the shapes 
and spins of the dark matter and gas components of the 10,000 most massive objects extracted from 
the simulation as well as the gas fraction in those objects. We find that the shapes of objects tend 
to be prolate both in the dark matter and gas. There is a clear dependence of shape on halo mass, 
the more massive ones being less spherical than the less massive objects. The gas distribution is 
nevertheless much more spherical than the dark matter, although the triaxiality parameters of gas 
and dark matter differ only by a few percent and it increases with cluster mass. The spin parameters 
of gas and dark matter can be well fitted by a lognormal distribution function. On average, the spin 
of gas is 1.4 larger than the spin of dark matter. We find a similar behavior for the spins at higher 
rcdshifts, with a slightly decrease of the spin ratios to 1.16 at z = 1. The cosmic normalized baryon 
fraction in the entire cluster sample ranges from YJ, = 0.94, at z = 1 to Yj, = 0.92 at z = 0. At both 
redshifts we find a slightly, but statistically significant decrease of with cluster mass. 
Subject headings: cosmology: theory - clusters: general - methods:numerical 



1. INTRODUCTION 

According to the standard scenario cosmic structures 
form by the gravitational collapse of density fluctuations. 
This collapse is mainly determined by the dark matter 
(DM) which contributes 85% to the total matter density 
in the universe. Baryonic matter follows the DM and 
forms visible objects like galaxies and clusters inside the 
DM halos. 

Shape and angular momentum are two important 
characteristics of halos. The shape of DM halos 
has been already widely studi ed, mainly by mean s 
of N-body si mulat io ns (e.g. lAllgood et~atl (|2006h . 
iMaccioet all (12006ft . iBett et all ?2006l) and references 
therein). I Avila- Reese et al.l ( 20051 ) discuss the depen- 
dence of the shape o f galaxy sized halos on environment. 
iBailin k, Steinmet3 (|2005( ) studied the internal shape of 
DM halos. Clusters of galaxies are the most recently 
formed objects in the universe. Most of the gas has not 
had time to cool. Since gas contributes only 15 % of the 
total mass one would expect that the gas follows DM and 
both distribution should be similar. In order to quantify 
this similarity, we have characterized the shape of both 
components in halos obtained from a large non-radiative 
cosmological gasdynamical simulation. Further we have 
compared the baryon fraction in clusters with the cosmic 
baryon fraction. 

The origin, evolution and distribution of DM halo 
spins have been wide l y disc us sed in the p a st als o 
(eg. IVitvitska et all (120021) iBullock et all (|2001h . 
iMaccio et all (|2006ft . iBett et all (|2006l )). In galactic size 
halos the angular momentum of the gas component is im- 
portant for the understanding of disk formation. Using 
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gasdynamical simulation s Ivan den Bosch et al.l (120021), 
IChen, Jing. &: Yoshikawl (|2003h . ISharma fc Steinmetd 
(2005) have studied the spin of the gas component in 
galactic halos and found that at redshift z — the gas 
component, on average, has a larger spin than DM. Here 
we will discuss briefly the relation between the different 
definitions of the spin parameter and study the distribu- 
tion of spin parameters of DM and gas components in 
our numerical clusters. 

Most of the previous studies on shape and spin of ha- 
los used either collisionless dark matter simulations with 
large number of particles and relatively large volumes or 
gasdynamical (dissipative and non dissipative) simula- 
tions in small volumes with a small number of objects. 
The advantage of our study is that we have both gas and 
dark matter components in the MareNostrum SPH sim- 
ulation and a large cluster sample (more than 10,000) 
to perform statistics, thanks to the large computational 
volume and the high number of particles of this simula- 
tion. 

2. NUMERICAL SIMULATION 

This simulation, dubbed The MareNostrum Uni- 
verse, was performed with the entropy conserving Gad- 
GET2 code (|Springelll2005l ) on the MareNostrum super- 
computer at the Barcelona Supercomputer Center using 
the equivalent of about 29 years on a single CPU. It fol- 
lowed the non linear evolution of structures in gas and 
dark matter (DM) from z — 40 to the present epoch 
(z = 0) within a comoving cube of 500/i _1 Mpc edges. 
We assumed the spatially flat concordance cosmological 
model with the following parameters: the total matter 
density fi m = 0.3, the baryon density = 0.045, the 
cosmological constant ilj^ = 0.7, the Hubble parameter 
h = 0.7, the slope of the initial power spectrum n = 1 
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Fig. 1. — Left: Shape of the dark matter distribution in clusters with total masses larger than 5 X 10 li h 1 Mq, a\, ai, az are the three 
main axes of the clusters. Right: The same for the shape of the gas distribution. 



and the normalization erg = 0.9. The power spectrum 
used to generate the initial conditions for the simulation 
was kindly provided by Wayne Hu in a tabulated form. 
It was obtained from a direct integration of the Boltz- 
mann code for the parameters described above. We did 
numerical interpolation to compute the contribution of 
the different Fourier modes. Both components, the gas 
and the DM, were resolved by 1024 3 particles, which re- 
sulted in a mass of 8.3 x 10 9 h~ 1 M Q for the DM particles 
and 1.5 x 10 9 /i _1 M Q for the gas particles, respectively. 

Gadget2 uses the TREEPM algorithm on a homo- 
geneous Eulerian grid to compute large scale forces by 
the Particle-Mesh algorithm. In this simulation we em- 
ployed 1024 3 mesh points to compute the density field 
from particle positions and FFT to derive gravitational 
forces. Within Gadget2 the equations of gas dynam- 
ics are solved by means of the Smoothed Particle Hy- 
drodynamics (SPH) method in its entropy conservation 
scheme. We did not include dissipative or radiative pro- 
cesses or star formation. The spatial force resolution was 
set to an equivalent Plummer gravitational softening of 
15 h^ 1 comoving kpc. The SPH smoothing length was 
set to the distance to the 40 th nearest neighbor of each 
SPH particle. 

The MareNostrum universe is part of a series of 
simulations we have performed during the last years. To 
this end, we had generated a realization of the ACDM 
power spectrum with 2048 3 particles and then decreased 
the mass resolution in the whole box to 256 3 , 512 3 and 
1024 3 particles respectively. The lower resolution (2 x 
256 3 , 2 x 512 3 ) SPH simulatio ns have been used to study 
properties of galaxy clusters (I Yepes et al.ll2004j) and the 
shape-alignment relation o f clusters (|Faltenbacher et all 
120051: iBasilakos et al.|[2Q06f ). The MareNostrum Universe 
have been recently used to analyzed the entropy profiles 
of the gas and DM in galaxy clusters (Faltenb acher et ail 
2006). Due to computational limitations, we could not 
yet simulate the evolution of the full box with the maxi- 
mal possible resolution of the initial conditions, 2 x 2048 3 
particles. However, using the multi mass technique de- 
scribed in iKlypin et alj (|2001f) selected individual clus- 



ters have already bee n re-simulated at this resolution 
(jAscasibar et al.ll2006l ). 

It is a challenge to find within a distribution of 2 billion 
particles all structures and substructures and to deter- 
mine their properties. Here we have used a newly devel- 
oped parallel versi on of the hierarchic al friends-of-friends 
(FOF) algorithm (|Klypin et al.lll999f ). In a first step we 
construct the minimum spanning tree for the distribu- 
tion of gas and DM particles. After topological sorting 
we get a cluster-ordered sequence from which we can eas- 
ily extract FOF clusters by simple cutting the tree at the 
desired linking length. We use a basic linking length of 
0.17 of the me an interparticle separat ion to extract the 
FOF clusters (|Gottlober et all l2006bT) . We divide this 
linking length by 2™ (n = 1,3) to find substructures and 
in particular the centers (density peaks) of our objects. 
We were running the minimum spanning tree and the 
FOF analysis independently over DM and gas particles 
to find their distribution. Using a linking length of 0.17 
at redshift z = 0we have identified more than 2 million 
objects with more than 20 DM particl es which closely 
follow a Sheth-Tormen mass function ((Gottlobe r et al.l 
I2006al) . 

To determine the shape and spin of the clusters we 
have selected a subsample of more than 10,000 clusters 
and groups with masses larger than 5 x 10 13 /i _1 M©. The 
lower mass threshold corresponds to clusters rsp. mas- 
sive groups with about 5000 gas particles and 5000 dark 
matter particles. As mentioned above the objects were 
identified independently from the gas and DM distribu- 
tions. Due to the large number of particles in this sim- 
ulation we can unambiguously match DM and gas halos 
with the same center of mass to one cluster. In the fol- 
lowing we will study these two components of the clusters 
in more detail. 

3. SHAPES 

Using the FOF method one extracts rather complex 
objects from the simulation which are characterized by 
an iso-density surface given by the linking length. In 
first approximation these objects can be characterized 
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Fig. 2. — Ratio of the minor to major axes a^/ai (filled symbols) 
and triaxiality parameter Z (open symbols) for the DM (circles) 
and gas (triangles) distribution in our numerical clusters 

by three-axial ellipsoids. The shape and orientation of 
the ellipsoids can be directly calculated as eigenvectors 
of the inertia tensor of the given object. Then the shape 
is characterized by the ratios between the lengths of the 
axes ai > ai > a 3 . 

We ran the FOF algorithm independently over the DM 
and gas particles to determine the shape of the distribu- 
tion of the two components. In Fig. [1] we show the 
shape of the dark matter and gas distribution in clusters 
with total masses larger than 5 x 1O 13 /i _1 M0. The ra- 
tios 03/02 and a-i/a\ of the clusters have been sampled 
into 25 bins of size 0.04. The plot shows a clear difference 
between the shapes of DM (left) and gas (right) distribu- 
tion. DM halos are centered at ratios (0.7, 0.75) whereas 
for the gas halos the center is at (0.85, 0.9), {i.e. the 
gas halos are much more spherical). Both shapes tend to 
be more prolate, i.e. 03/02 > aija\. This could be the 
result of mergi ng along a preferred dire ction, the large 
scale filaments (jFaltenbacher et aLll2005j ). 

The exact position of the maximum in left and right 
panel of Fig. Q] depends on the low mass cut-off in the 
cluster sample, but the qualitative behavior is indepen- 
dent of it. In [21 we show how the ratio of the minor to 
major axis 03/01 of DM (filled circles) and gas (filled tri- 
angles) halos depend on cluster mass. Note, that the dis- 
persion in the values of the ratios 03/01 is large (within 
one bin the ler scatter is 0.14). This can be seen in 
Fig. [1] for the ratios 02/01 and 03/02- However, due 
to the large number of objects per mass bin the stan- 
dard deviation of the average within each bin shown in 
Fig. [5] is small and sometimes smaller than the symbol 
size. The solid line is a linear fit in the semi-logarithmic 
plot 03/01 = b + cln(M), where the mass M is given 
in h~ 1 Mo and b and c are (1.904, -0.0426) and (1.932, 
-0.0401) for DM and gas, respectively. A similar trend 
(but h igher 03/01 have been found bv iKasun fc Evrardl 
(|2005| ) in the Hubble Volume N-body simulation. 

The cluster shape can be further characterized by the 
tri axiality p aram eter Z = (ai — 02)/ (ai — 03) introduced 
bv iBinnevI (|1985h . Z = 1 corresponds to prolate ellip- 
soids of revolution (prolate spheroids), while Z = cor- 
responds to oblate ones. In Fig. [2] the triaxiality param- 



eter Z is shown for the DM (open circles) and gas halos 
(open triangles). The dotted line is a fit Z = d + f\n(M) 
with (d, /) = (-0.61, 0.039) for DM and (-0.43, 0.034) for 
gas halos. Contrary to the ratio 03/01 there is only lit- 
tle difference between the triaxiality parameter Z of the 
DM and gas halos. The younger the halos are (i.e. more 
massive) the more prolate they look. 

4. SPIN 

In N-body numerical simulations the identified clusters 
can be characterized by their mass velocity and angular 
momentum or spin. O r iginal ly the spin parameter A was 
introduced bv iPeeblei (|1971h as the "convenient dimen- 
sionless number" 



A = 



JjE] 1 / 2 
GM b / 2 ' 



(1) 



where J is the total angular momentum of the object, E 
its tot al energy and M its mass. Following lPadmanabhanl 
(1993) one can interprete A as the ratio of the angular ve- 
locity oj of the system to the angular velocity w S up of the 
system that would provide rotational support. Charac- 
terizing the system by its angular momentum J ~ luAIR 2 
and rotational support by ui 2 up R 2 



GM/R one finds 



A 



J 



-'sup 



G l/2 M 3/ 2jR l/2' 



(2) 



which is up to a factor of y/2 Bullock's A' (|Bullock et al.1 
l2001h if one replaces GM/R by the circular velocity V 2 irc . 
Contrary to Eq.fll]) both Bullock's definition and Eq.([2|) 
contain only quantities which can be easily calculated in 
numerical simulations for spherical halos. If we assume 
that the total energy of the system is characterized by 
E ~ —GM 2 /R one comes back to Peebles' original defi- 
nition. 

The halos in the numerical simulation are assumed to 
be virialised and thus characterized by 2T+U = 0, where 
T an d U are the kinetic and potential energies of the ob- 
ject. iBett et al ()2006l ) show in their Fig. 5 how real 
halos scatter around this assumption. Replacing the po- 
tential energy in Eq.JT]) by (— 2T) we end up with a third 
definition of the spin parameter: 



A 



JT 1 ' 2 
GM 5 / 2 ' 



(3) 



As with the previous definition this A can be easily calcu- 
lated numerically. The advantage of this spin parameter 
is that it can be calculated easily for any shape of the 
halo. Therefore, this definition is especially suited for 
the calculation of the spin parameter of halos found with 
the FOF analysis. 

The halos identified in the MareNostrum universe con- 
sists of two components, dark matter and gas. Then for 
each component the spin parameter is the ratio of the an- 
gular velocity u> of that component to u sup of the system. 
Therefore, we have 



A 



gas(DM) 



-gas(DM) ,.s 

s(DM) ^(Afgas+MDM^vir) 1 / 2 ' ' ' 

denote the gas (DM) mass inside the 
virial sphere of radius i? V ir and J gas (DM) the angular mo- 
mentum of the corresponding component. Note, that we 
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Fig. 3. — Comparison of the DM spin parameter A3 calculated 
according to Eq. (J3j and A4 computed from Eq. (J4j 

have introduced in Eq. (|4]) an addition al quotient of y/2 t o 
b e consistent with the d efiniti on of [Bullock ct al. ( 20011). 
Ivan den Bosch et al.l (|2002h and ISharma "fc" Stcimnctz 
discussed in detail the angular momentum distri- 
butions of the gas and dark matter components of halos 
found in SPH cosmological simulations. To this end they 
decomposed the particle velocities into a streaming and 
a random component. The velocities of the SPH gas 
particles are streaming velocities to which one could add 
microscopic random motions. The velocities of the DM 
particles contain random motions which could be elimi- 
nated by smoothing over a certain interpolation kernel. 
Nevertheless, as long as the number of particles in the 
halo is very large compared to the number of smoothing 
neighbors the effects of the smoothing on the spin param- 
eter can b e neglected (Sharma: private communication) . 
Following Shar ma &; Steinmeti (|2005f ) we use the total 
velocities of the DM and gas particles as provided by 
the simulation to calculate the spin parameters of both 
components. 

We want to compare the spin parameters of a given 
structure calculated by the different abovementioned def- 
initions. The objects found by the FOF algorithm have 
arbitrary shapes and the spin parameter as defined in Eq. 

A3 is calculated with respect to the center of mass of 
each object. On the other hand, we have also calculated 
the spin parameter for spherical halos at virial overden- 
sity (334 x mean density at redshift z = 0). In this case, 
we first identify for each of the FOF objects the highest 
density peak. To this end, we decrease the linking length 
by a factor of 8 ending up with sub-structures of about 
170.000x mean density. For each FOF object we take 
the center of mass of the most massive sub-structure as 
the center of a sphere and determine the virial radius 
and mass within this sphere. Then we use Eq. ^ to 
calculate the DM spin parameter A4. In Fig. [3] we com- 
pare the DM spin parameters A3 and A4. Having in mind 
the completely different treatments of the spin parame- 
ter they agree surprisingly well. More than 63 % of our 
halos are within 20 % scatter around the A3 = A4 rela- 
tion. The reason for the larger deviations shown in this 
plot is due to the fact that halos defined from the spher- 
ical overdensity criterium and those defined from FOF 



algorithm look completely different if substructures are 
present. In such cases, the spin of the FOF objects is 
calculated with respect to the center of mass whereas 
for the spherical halos it is calculated with respect to the 
the position of the highest density peak. These positions, 
and the corresponding spins may differ substantially. For 
those halos in which this effect is not important, the spin 
parameters defined by both methods have very similar 
values. However, there is also another source of scatter, 
although much smaller than the previous one, due to the 
assumption of virial equilibrium IT + U = which is 
not exactly fulfilled in our halos. Thus, even for spheri- 
cal o bjects the spins A3 a nd A4 will not be identical. In 
fact. lShapiro et ail (|2004D pointed out that the presence 
of infalling matter acts as a surface pressure even for col- 
lisionless DM. The surface term leads to 2T + U > 0. 
The spin parameter calculated from the kinetic energy, 
A3 tends to be slightly larger than the spin parameter A4, 
probably because the total kinetic energy of DM parti- 
cles is larger than that assu med by virial equilibrium. 
iHetznecker fc Burkertj (|2006| ) have studied the distribu- 
tion of the virial coefficient 77 = 2T/\U\ for a sample of 
dark matter halos from ACDM N-body simulations and 
showed that the mean value of 77 approaches unity at 
redshift z = 0. 

In § [3] we found different shapes (and therefore vol- 
umes) for the DM and gas components of the halos. 
Therefore, to compare the spin parameter of both com- 
ponents we better used A defined by Eq. 2] We put 
a sphere at the position of the most massive substruc- 
ture as described above and determine the virial radius 
and mass of the halo as well as the gas and dark matter 
masses (see §[S]) and the corresponding angular momenta. 

The resulting distribution of the spin parameter of the 
gas and dark matter components are shown in Fig. 3J 
The distribution of the spin parameter A can be described 
by a log-normal distribution 

™ = ^-ex P r ln2(A/A ° r 



/2na x X ' - (5) 

The best fit parameters are for the DM distribution 
A = 0.0351±0.0016,cr A = 0.6470±0.0067andforthegas 
distribution A = 0.0462 ± 0.0012, a x = 0.6086 ± 0.0030. 
P(A) has a maximum at X max — 0.02 31, 0.0319 f or the 
DM resp. gas distribution. Recently, iBett et al (|2006l ) 
have proposed another distribution function which fits 
their data better. Since the scatter of Ao with varying 
total number of halos depending on the lower mass cut- 
off is in our case of the same order as the error in the 
determination of Ao we believe that the log-normal func- 
tion is a sufficient fit to describe the behavior of the spin 
parameter. The differences of Ao using different defini- 
tions of A are at least of the same order. 

Since dark matter dominates the total mass the total 
distribution of the spin parameter of dark matter and 
gas practically coincides with that of the dark matter 
only. Moreover, it is in agreemen t with the distribution o f 
the spin parameter calculated bv lGottlober et all ((2006a) 
u sing Eq.flSl) for the non-sphe rical FOF halos. 

Ivan den Bosch et al.l (|2002f ) found agreement between 
the spin distribution of the gas and DM components 
of 378 halos identified in a small box of 10/i~ 1 Mpc 
at redshift z — 3. We see a substantial shift of the 
log-normal distribution of the gas spin towards higher 
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Fig. 5. — Distribution of the spin parameter of the gas and DM 
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Fig. 4. — Left: Distribution of the spin parameter of the dark matter in halos with masses larger than 5 X 10 13 /i _1 Mq. The dotted line 
is a log-normal distribution Eq.JS} with the parameters Ao and a as given in the plot. Right: The same for the spin parameter of the gas 
in halos. 

mine the baryon fraction is by sure a simplification, but 
since we are interested in the baryon fraction at virial ra- 
dius we expect not to be affected very much by neglecting 
cooling and star formation processes (see the discussion 
below) . 

During the last decade baryon fractions in clusters 

Yfr /cluster //cosmic /cluster/ (^fe/^matter) have been 

studied in non-radiative simu lation by several authors. 
lEke. Navarro, fc Fr cnk (1998}) found within the virial ra- 
dius a baryon fraction of 0.85 - 0.9. Within the Santa 
Barbara cluster comparison project (|Frenk et al.l Il999h 
the baryon fraction at virial radius averaged over all 
codes was Y& = 0.92. They found a substantial scatter 
between codes and a systematic offset between SPH and 
grid codes which led to higher b aryon fractions. Recently, 
iKravtsov. Nagai. fc Vikhlininl (|2005f ) found that at large 
radii the baryon fraction in the ART simulations is by 
about 3% - 5% higher than in GADGET simulations of 
the same clusters. Within our non-radiative simulation 
we have explored the baryon fraction at the virial radius 
in objects with virial masses larger than 5 x IO^/i^Mq. 
The assumed virial overdensity at redshift z = is 334, at 
redshift z = 1 it is 201. At redshifts z — we found more 
than 10000 objects and about 2500 at redshift z = 1. The 
scatter in the measured baryon fraction of our cluster was 
found to be quite large. At redshift z = it ranges be- 
tween 0.85 and 1.0 with a mean of 0.92 (fGottlobe r et al.l 
l2006ah . 

In Fig. [S]we show the ratio Yf, as a function of the virial 
mass of clusters. Due to the large number of objects the 
standard deviation of the mean value of Yj, for different 
mass bins is small and we could fit a linear relation in 
the semi-logarithmic plot, Yj, = aln(M) + (3, with the 
slope a — —0.005 ± 0.001 at redshift z — and a = 
-0.006 ±0.001 at redshift z = 1. There is a 2% decrease 
of the baryon frac tion between redshifts one and zero. 
lEttori et al.1 (j2006) found 5 % decrease of the baryon 
fraction in this time interval. Within our model the de- 
crease can be naturally expected. Since the gas works 
against its pressure during the formation of the cluster 
the total dark matter mass grows slightly faster than the 
gas mass. Therefore, the relative gas fraction decreases 
with time as can be seen in Fig. [6] For illustration 
we have selected 1200 clusters with masses larger than 
2 x 1O 14 /i _1 M0 and 2900 clusters with masses smaller 
than 2 x 10 14 ft.~ 1 M Q and compared the mass growth of 



spin in comparison to the dark matter. Our Ao are 
sli ghtly larger than the values (0.0287, 0.0412) reported 
by ISharma fc Steinmetzl (|2005l ) for a sample of 41 ha- 
los in a box of 32.5/i _1 Mpc at redshift z = 0. We 
find a mean \ ga s/^DM of 1.39 (see Fig. [5]) with a 
tendency to decrease with halo mass. The standard 
deviation is 0.57. Only the 159 most massive halos 
(M vir > 5 x 10 14 ft. _1 M Q ) show a substantially smaller 
^gas/^DM = 1-23 with a standard deviation of 0.45. 

At z = 1 we found a similar behavior of the spins 
of the DM and gas components. Both spins follow log- 
normal distributions which can be described by Ao — 
0.0467 and 0.0541 respectively, i.e. the mean spin of 
the gas component is 1.16 times larger than the spin of 
the DM component. During further evolution this ratio 
increases until 1.39. The spins themselves are slightly 
larger at redshift z = 1. 

5. BARYON FRACTION 

The observed cluster baryon fraction is an important 
tool for the determination of cosmological parameters. 
Typically the gas fraction in clusters is measured by X- 
ray observations at overdensities 500 or larger, i.e. well 
inside the virial radius. A certain fraction of the baryons 
reside in stars. Using a non-radiative simulation to deter- 
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Fig. 6. — Baryon fraction, normalizwed to the cosmic baryon 
fraction, in clusters as a function of halo mass for two different 
redshifts. 

their DM and baryonic components between z = and 

Z = 1, D DM, gas = m D M,gas{z = 0) / m DM, gas (z = 1). 

The averaged relative growth DdmI 'D gas is 1.03 for 
M > 2 x 10 u /i- x M Q and 1.02 for M < 2 x 1O 14 /i- 1 M 0) 
i.e. clusters accrete DM faster than gas and the mean 
baryonic fraction decreases with time. Originally the 
baryon fraction was homogenous. During collapse of 
the whole cluster region the baryonic fraction decreases 
slowly. Most massive clusters are situated in the posi- 
tions of the highest density peaks in the originally al- 
most homogenous medium. These peaks start to grow 
first. This could be the reason for small but statisti- 
cally signific ant slope a which w e found. Note, however, 
that recently lCrain et al.l (|2006h have found the opposite, 
a slightly increasing baryon fraction with mass, from a 
set of clusters extracted from a non-radiative SPH re- 
simulation of the Millenium Run, but with 2 times less 
number of particles than ours. Also, the relative baryon 
fraction of their clusters (at overdensity 200) falls below 
0.9. 

Recently iKravtsov. Nagai. fc Vikhlininl (|2005D studied 
the effects of cooling and star formation on the the 
baryon fractions in clusters. They showed that the cool- 
ing of the gas and the associated star formation increases 
the baryon fraction even at radii as large as the virial 
radius. The averaged baryon fraction at virial radius 
differ by 5 % between the simulations with cooling and 
star formation (Yb = 1.02) and the non-radiative simu- 
lations (Yb = 0.97), which is already ~ 5 % larger than 
the value we obt ained from the Mar eNostrum Universe 
SPH simulation. lEttori et al.l (|2006f ) have also studied 
a set of clusters obtained both in non-radiative SPH 
simulations as well as in simulations including cooling 
and star formation. For the non-radiative simulations 
they found an averaged baryon fraction of 0.89, which is 
slightly lower than the value w e found . In agreement with 
IKravtsov. Nagai. fc Vik hlinin (2005) they also found in 



their simulations with cooling and star formation that 
the baryon fraction at virial radius increases with respect 
to the non-radiative case by about 3 %. 

6. CONCLUSIONS 

We have selected a statistically significant sample of 
more than 10,000 clusters from the MareNostrum uni- 
verse simulation. With this database we have determined 
the shape and the spin of the DM and gas components of 
the clusters independently. We found that both the gas 
and dark matter components tend to be prolate although 
the gas is much more spherically distributed. The mean 
ratio of the minor to major axis decreases with increasing 
halo mass, i.e. younger objects tend to be less spherical , 
in agreement with the results of iKasun fc Evrardl (2005) 
obtained from the Hubble Volume DM only simulation. 
On the other hand, the triaxiality parameters of the gas 
and DM component differ by a few percent only and both 
of them increase with halo mass. 

The spin parameters of the DM and gas components 
are well represented by lognormal distribution functions 
with A = 0.0351 (DM) and A = 0.0462 (gas). These 
values are slightly larger than the values found for a sam- 
ple of 41 halos with masses from dwarf to bright galaxies 
(jSharma fc Steinmet j |2005). On average, the spin of the 
gas component is by a factor of 1.39 larger than the spin 
of the DM component. This ratio is smaller for very mas- 
sive halos, 1.23 for halos with M vir > 5 x 10 14 /i _1 M o . 
It decreases slightly with redshift being 1.16 at z = 1. 

On average, the baryon fraction in the cluster sample 
is Yb — 0.92. The baryon fraction increases with redshift, 
being 0.94 at z = 1. At both redshift the baryon frac- 
tion decreases slightly with increasing virial mass of the 
clusters. The radiative processes of cooling and star for- 
mation are expected to change the baryon content found 
for halos in non-radiative simulations. There is a discrep- 
ancy between the predicted baryon fraction in SPH and 
AMR codes. Both effects are of the same order as the 
time and mass dependence found here. Thus, at present 
observational projects to estimate cosmological parame- 
ter such as the total matter density in the universe and 
the equation of state of the dark energy can still safely 
assume a rather universal, unevolving baryon fraction for 
clusters of galaxies. 
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